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MHD Dynamical Relaxation of Coronal Magnetic Fields 

I. Parallel Untwisted Magnetic Fields in 2D 
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ABSTRACT 

Context. For the last thirty years, most of the studies on the relaxation of stressed magnetic fields in the solar environment have only 
considered the Lorentz force, neglecting plasma contributions, and therefore, limiting every equilibrium to that of a force-free field. 
Aims. Here we begin a study of the non-resistive evolution of finite beta plasmas and their relaxation to magnetohydrostatic states, 
where magnetic forces are balanced by plasma-pressure gradients, by using a simple 2D scenario involving a hydromagnetic dis- 
turbance to a uniform magnetic field. The final equilibrium state is predicted as a function of the initial disturbances, with aims to 
demonstrate what happens to the plasma during the relaxation process and to see what effects it has on the final equilibrium state. 
Methods. A set of numerical experiments are run using a full MHD code, with the relaxation driven by magnetoacoustic waves 
damped by viscous effects. The numerical results are compared with analytical calculations made within the linear regime, in which 
the whole process must remain adiabatic. Particular attention is paid to the thermodynamic behaviour of the plasma during the relax- 
ation. 

Results. The analytical predictions for the final non force-free equilibrium depend only on the initial perturbations and the total pres- 
sure of the system. It is found that these predictions hold surprisingly well even for amplitudes of the perturbation far outside the 
linear regime. 

Conclusions. Including the effects of a finite plasma beta in relaxation experiments leads to significant differences from the force-free 
case. 
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1. Introduction 

The magnetic field of the solar cor ona is believed to evolv e 
through a series of force-free states dHevvaerts & Pries3ll984l) . 
Since the solar corona involves a low-beta plasma in which 
magnetic forces dominate over plasma forces, this is not 
an unreasonable assumption, and so, most of the recent 
studies on th e relaxation of co r onal magnetic fields (e.g. 
Mackav & van Ballegooiie nI l2006t iJi et al. l 120071: llnoue et al.1 
20081; iJanse & Lowii2009i: Miller et al.ll2'009t iPontin et al.ll2009l) 



have been done by considering the approximation of an ex- 
tremely tenuous plasma, for which the plasma pressure does not 
play an important role, and the persistent hydromagnetic struc- 
tures of the solar corona are assumed to be in magnetic balance, 
with zero pressure gradients. On the other hand, only within the 
past few years, Ruan et al. (2008) and Gary (2009) have started 
to consider the reconstruction of the global coronal magnetic 
field including a finite Lorentz force balanced by magnetic and 
gravity forces. 

In addition, there are many codes available to calculate 
those force-free fields from the observed magnetic field in the 
photosp here ( Amari et al. 1998; Wiegelmann et al. 2006, 2008; 
ISchriiver et all 120061: iMetcalf et al.i [2008). These codes have 
been used with varying degrees of success to determine the mag- 
netic fi^ld_of^olarflaresa De Rosa et al. 
2009; 'Sc hriiver et al.°2008VReg nied2008HWhea tland & Regnierl 
,2009). However, problems remain with these approaches. In par- 
ticular, a non-linear force-free field determined from a line-of- 
sight photospheric magnetic field is not unique, but is one of an 



infinite number of possible solutions. This fa ct is well k nown 
and has been discussed by several authors (see lLowll2006l) . 

Also, the low beta plasma assumption is only valid for the 
solar corona, but is not valid in the photosphere, chromosphere 
or much of the transition region where plasma effects become 
more important. Furthermore, in the solar corona, the tenuous 
plasma will be able to create a high beta in the surr oundings of 
a magnetic null point (iMcLaughlin & Hoodll2006l) . Moreover, 
even where the plasma has a low beta, there are still some ef- 
fects on the final equilibrium of the magnetic field that will lead 
to energetic consequences as the field is relaxing to its new equi- 
librium state. 

In particular, relaxation can involve the magnetic field evolv- 
ing from a stressed state to a force-free field with a lower 
energy, and for a consistent scenario where the energy can- 
not escape the system, conservation of energy implies that the 
loss es of magnetic e nergy must be conver ted into something 
else. lBrowning et a T (2008) and Hood et aH (1200 9) investigated 
Taylor relaxation (jTaylor .1974.) through a series of non-linear 
3D simulations, initiated by an ideal MHD instability. Although 
the initial state was force-free, the final state involved a high 
temperature plasma with a significant value for the plasma beta. 
Placing aside some extra contributions such as radiative losses, 
if a substantial fraction of the magnetic energy released goes into 
the internal energy, then the plasma beta cannot be small. Hence, 
considering the behaviour of the plasma is important even if it 
has li ttle eff ect on the final magnetic equilibrium. On the other 
hand. iGarvl (j2001) suggested the possibility that there is high 
beta plasma in the solar corona above active regions. 
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In this paper, we will start by considering a simple scenario, 
with a uniform 2D magnetic field. There will be no current sheet 
formation and no changes in the magnetic topology. The aim 
of the paper is to investigate the relaxation of the hydromag- 
netic fluid for various diff'erent values of the plasma beta. To 
perturb the system, we introduce a local small enhancement (or 
deficit) in the plasma pressure (or in the magnetic field), un- 
der the frozen-in condition. The subsequent relaxation leads to a 
new magnetohydrodynamic equilibrium in which the perturba- 
tion has been dissipated by viscous forces acting on the flow. 

It is worthwhile mentioning that relaxation via Ohmic dis- 
sipation, due to the effect of resistivity, or magnetic diffusiv- 
ity, represents a substantially different problem; While viscos- 
ity dissipates the plasma velocity, diffusivity tends to eliminate 
the electric current density, and such a relaxed state can only 
involve potential fields, which are mathematically well defined 
and are uniquely determined by the components of the magnetic 
field normal to the boundaries. Furthermore, the time-scales for 
an Ohmic relaxation in very high magnetic Reynolds number 
environments, such as the solar corona, are in general probably 
larger than the age of t he S un itself, outwith regions with very 
small length scales (see lPriesllll982h . 

By considering the effects of the plasma pressure in the re- 
laxation, we are facing a totally different problem from that of 
the force-free relaxation studied by many others. In our non- 
resistive MHD relaxation, the plasma displacements driven by 
the initial pressure enhancement will carry the magnetic field 
with them, generating an electric current and a magnetic pres- 
sure. Hence, the resulting equilibrium will have to involve a bal- 
ance between the Lorentz force and the plasma-pressure gradi- 
ent. 

The effects of including a finite plasma beta are relevant not 
only in the high plasma beta regions of the solar atmosphere such 
as the photosphere and chromosphere, but will also be relevant 
in the solar corona. Obvious regions where the plasma beta is 
likely to have a significant effect are in the vicinities of mag- 
netic neutral points, where the magnetic field vanishes. These 
configurations will be the subject of a further paper, while in 
the present paper we study the plasma beta effects on the sim- 
plest magnetic configuration, a uniform magnetic field, which is 
absolutely general and might be compared with different solar 
environments such as a region in a coronal prominence or part 
of a coronal loop. 

In Section |2l we first present the setup of the two- 
dimensional linear problem. In Sections 12.11 and 12.2! we solve 
the equations for a one-dimensional vertical and horizontal per- 
turbation, respectively, and in l2.3l we obtain the whole solution 
from a general two-dimensional linear perturbation, producing 
a practical and precise analytical solution to the problem. In 
Section [3] we presen t a set of numeri cal experiments run us- 
ing the Lare2D code (' Arber et"al]|2001h . and we compare these 
with the previous analytical results. The effects of the non-linear 
terms are considered when the magnitude of the perturbation is 
increased, in Section|4l followed by the summary and some con- 
clusions which are presented in Section|5] 



2. Linear 2D Equations 

The initial set up involves a uniform magnetic field pointing in 
the vertical y-direction, Bq = BoC, , and a background plasma 
of a constant gas pressure, density and temperature. The ini- 
tial disturbances are supposed to be small, in order to stay in 
the linear regime. Expressing each quantity q{x,y, t) as the sum 
of a background constant value plus a perturbation, q{x, y, t) - 



qo + qiix,y,t), and substituting those expressions into the vis- 
cous, non-resistive, two-dimensional MHD equations, with no 
gravity force for simplicity, and neglecting terms involving prod- 
ucts of perturbations, we obtain the following set of linearized 
equations, in scalar form: 



Po 



Po 



dp\ 
dt 
dv^ 

dt 

dpi 
dt 

dB, 
dt 

dB, 



-poV • vi , 

dpi Bo dBy Bo dB, 



(1) 



. . +— ^+Pov Vn',+ -— (V-v.) ,(2) 

ox pq ox Po dy \ i ox 



1 d 



dpi i„2 1 5 , 

-ypoV ■ Vi , 

„ dvx 



-B, 



dy ' 

dvx 
dx 



(3) 
(4) 
(5) 
(6) 



where po, po and Bq are the constant background plasma den- 
sity, plasma pressure and magnetic field, v is the kinematic vis- 
cosity. Pi, Pi and Vi are the perturbations on plasma density, 
plasma pressure and velocity, and v^, v,,, B ^ and By are the x and 
y components of the perturbed velocity and perturbed magnetic 
field, respectively. Also, the plasma pressure, p, density, p, and 
internal energy, e, are related by the perfect gas law, for an ideal 
polytropic gas. 



p = pe(y - 1) 



(7) 



where y is the ratio of specific heats, often assumed to be 5/3 for 
a highly ionised hydrogen plasma. 

In a general scheme, viscosity would add a heating term to 
the energy equation, but this term is second order, and so, the 
process is adiabatic within the linear regime, and there is no heat- 
ing of any kind taking place. Hence, the entropy per unit mass, 
S = p/p^, is conserved, for each single fluid element, and for 
the entire box. 

From the conservation of entropy, a relation between the 
plasma pressure and density perturbations may be obtained, 
within first order (i.e. neglecting the terms involving products 
of perturbations): 



Po + Pi 
Hence, 



PO L ^ Pi ypi 

Po I po po 



constant 



Aj7l 
PO 



po 



where A indicates the difference between final and initial state of 
the perturbation, such that 



Api 



c?Api 



(8) 



where Cj = yjjpolpo is the sound speed. 

To solve these equations, we first determine the solution for 
a one-dimensional perturbation, which depends only on x (Sec. 
12.11 ). then only on y (Sec. 12.2b . and finally, using the ID results, 
we derive the solution for a more general 2D perturbation (Sec. 
|231). 
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2.1. Perturbation dependent on X 

Consider a perturbation varying only in the direction perpendic- 
ular to the magnetic field lines, x. The magnetic field vector will 
have a non-zero y-component, Bi(x, f) - By{x,t)ty, while the 
velocity will have a non-zero jc-component, y\{x, f) - Vx{x, f)e,i. 
Equations ([T]i to (|6l) reduce to 



dpi 
dt 


dv^ 
ox 


(9) 




dpT yd^v^ 
dx dx^ 


(10) 


dpi 
dt 


dv, 

= -7Po-^ , 

ox 


(11) 


dBy 

~dr 


B 

dx 


(12) 



with v' - 4v/3, where v is the kinematic viscosity, and pr is the 
perturbed total pressure, given by 

BoBy 



Pt = Pi + 



(13) 



The equation governing the final equilibrium state can be ob- 
tained using Eq. ( fTOb . At equilibrium, the time dependence dis- 
appears, and the velocity is zero, thus, the equilibrium must have 
constant total pressure; 



dpT 
dx 



= . 



(14) 



Combining Equations ( fTTT i and ( fT2l l. we get the evolution of the 
total pressure as 

1 dpT 
Po dt 



-{cj + cl) 



2^dVx 

dx 



where Cs is the sound speed, defined above, and ca - BqI ^J/iqPo 
is the Afven speed. By the appropiate combination of Equations 
( [Tol l and ( fTsT i we get the wave equation for fast magneto-acoustic 
waves for the total pressure: 

d^PT 



dt^ 



(c. + c\) 



2^d^PT , , d I d^pT 



dx^ dt 



dx^ 



(16) 



Assuming that the total pressure can be considered as a contin- 
uous, periodic function, the solution of the last equation can be 
expressed as a superposition of plane waves, such as 



Pt(x, t) = Ke 



y k 



<Pke 



i(kx-iAjt) 



(17) 



where each wave number k corresponds to a different oscillation 
mode and is associated with a complex frequency a>(k) given by 
the dispersion relation, 

a/ + ikWcj - (c^ + c\)k^ = . 

The dispersion relation has the solution cli - a - bi, where a is 
the real frequency of the wave, and b is the damping term: 



■ c\) - k^v'^ 



b = -k^v' 
2 



In order to have a harmonic mode, the wave number k must sat- 
isfy k^V^ < 4(Cj + c^), and higher modes will be damped with- 
out any type of oscillation. On the other hand, notice that w = 



when k - 0. The undamped mode k - corresponds to the con- 
stant Fourier coefficient in the expansion of Eq. ( fTTb . which does 
not change in time, and is given by Fourier analysis as the homo- 
geneous redistribution of the initial total pressure to its average 
value. As f — > oo, it is only this constant that remains, as all the 
other terms are proportional to e Hence, this homogeneous 
redistribution is exactly the constant perturbed total pressure that 
defines our final equilibrium state: 



PT(oa) 



B^Byix,0) \ ^ 
p 1 (X, 0) H I ax . 



(18) 



where is the length of the x-domain. 

From the solution of Eq. ( fTSI l, and Eq. (fTsT i, we obtain an 
expression for V[(x,t), which, after substitution into Equations 
(01, (HD and O, gives 



dpi 1 dpT 

'dT ~ cl + cl~dr 

dpi _ c] dpT 
dt 
dBy 

'dT ~ po(c?+c2) dt 



Bo 



dpi 



(19) 



(20) 



(21) 



Integrating now from f = to f = oo, we obtain the perturbed 
quantities for the final equilibrium state, as functions of the per- 
turbed total pressure, to be added to the background values: 



Pi{x) = pi(x,0) + 



1 



(15) P\{x) = /?i(x,0)H- 



{pTioo)-pT(x,Q)), (22) 



ipT{^)-pT{x,0)), (23) 



Byix) = By(X,Q) + 



Bo 



poic] + c\) 



{PT{^) - Pt{x,Q)) . (24) 



Equations (l22l i. (l23T l and (l24l l state that no matter how we 
set our initial disturbance, the final equilibrium distributions are 
completely determined by the initial and the final total pressures 
of the system, which are given by the solution of the wave equa- 
tion. Note, that also the adiabatic equation for the linear regime 
given in Eq. ([8]l is satisfied. 

2.2. Perturbation dependent on y 

In the same way as above, we can get the solution for a pertur- 
bation varying only along the field. This time, we are dealing 
with a purely non-magnetic evolution, that will lead to a homo- 
geneous redistribution of the plasma pressure all along the field 
lines. The equations to solve are the y-components of Equations 
dlJ, (O and (|4|. The equilibrium is now given by 



dpi 

dy 



= 



and the solution for the perturbed quantities is 



pi(oo) 



{y,0)dy. 



1 



Pl(y) = pi(3',0)H- -r(/?i(oo)-;,i(3;,0)) , 



where L, is the length of the y-domain. 



(25) 

(26) 
(27) 
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2.3. Two-dimensional perturbation 

In this section, we study the relaxation of a general individ- 
ual two-dimensional perturbation within the linear regime. Once 
again, setting the velocities to zero, we get the equations govern- 
ing the final 2D equilibrium: 



B()By 

1^0 



BpdB, 
f^Q dy 
dpi 
dy 



, 




(28) 
(29) 



Equation (l29l l tells us that the final plasma pressure cannot 
depend on y, so the solution for the pressure must remain one- 
dimensional, as before. On the other hand, Eq. ( l28T l does not have 
a direct interpretation, as both spatial derivatives are involved. 
The term p\ + BqB^Iiiq represents the perturbed total pressure 
from Sec. 12.1! and the new term BoB^/fio represents the mag- 
netic tension due to the curvature of the field lines, which was 
zero in the ID cases. As usual, when looking for a periodic solu- 
tion, Fourier analysing makes life much simpler: Expressing our 
variables as functions of e'f'^ where each pair (k, I) repre- 
sents one single mode of oscillation in the global time evolution. 
Equations (l28l l and ( l29l l can be rewriten as 



PI 



where: 



BQBy\ R^^B, 



1^0 

I Pi 



i. The mode k - Q, 1 = represents the unperturbed back- 
ground values. 

ii. For k i^Q, 1 = 0, the equation of the equilibrium is 



—ipi + BqB,) = . 
ox 



(30) 



These modes only depend on x, and represent the homoge 
neous redistribution of the total pressure studied in Sec. |2TT 
iii. For k = 0, Z 9^ we get 



dpi 
dy 
dB, 
'di 



= 
= 



(31) 
(32) 



These modes do not modify By, instead they simply remove 
both the vertical gradients of magnetic tension and plasma 
pressure as in Sec. 12.21 Each of them is treated individually, 
as they are not coupled in the equations, 
iv. Finally, for those modes with k 0, / 0, we get 

kBy - IB^ = , 

which can be combined with the solenoidal condition for 
the magnetic field, V ■ B = 0, or, within our fourier notation, 

kBy + IB^ = . 

From these equations, we can conclude that, in the final 
equilibrium, the existance of a variation of By in the x- 
direction is totally incompatible with a variation of B^ in 
the y-direction. Hence, the modes with both wave numbers 
k and / non-zero may appear in the dynamical evolution, but 
not in the final equilibrium distributions. 



Therefore, with our uniform background magnetic field 
pointing straight in the vertical y-direction, the final equilibrium 
state is a combination of the background values (fe = 0, I = 0), 
plus the vertical non-magnetic evolution to a state with plasma 
pressure that is constant along y, and/or the smoothing of the 
horizontal component of the magnetic field (k = 0), plus the 
one-dimensional hydromagnetic evolution across the field lines 
(/ = 0). Note, that the perturbed B^ in the vertical direction (i.e. 
curved magnetic field) is not coupled with either pi or By, so 
the final magnetic field remains as straight lines, and By is not 
involved in the evolution of the total pressure of the system. 

Hence, we calculate the analytical 2D equilibrium in two 
steps: Firstly, the non-magnetic evolution in the vertical direc- 
tion. 



pKx) 



(x, y, 0) dy . 



1 



pl(x,y) = pi{x,y,0)+ -r(pl{x)- pi(x,y,0)) , 



(33) 



(34) 



and secondly, the hydromagnetic evolution in the horizontal di- 
rection, across the field. 



1 

pi(x,y) = pi(x,y)+— j(pt(°°) - Pt(x)) , 



pi{x,y) = pI(x) + ' (pt(°°) - Pj(x)) 



By(x,y) = By(X,0) + 



B, 







, .-(Pt(°°)- Pt(x)) ■ 



P0(C? H- C^) 



(35) 
(36) 
(37) 



In Equations ( 1331 ) to ( l37l i the quantities with a superscript repre- 
sent the state after the vertical evolution, with p*j{x) = p\{x) + 
Bofi,.(x,0)M). 

Looking back at the equations, we see a well known result 
from magnetohydrostatics, namely: In equilibrium, and in the 
absence of gravity, the plasma pressure must be constant along 
the field lines. The constant plasma pressure in the y-direction 
given by Eq. (l29l l is aligned with the straight magnetic configu- 
ration, with no magnetic tension, for the final equilibrium state. 

When analyzing the validity of the results above for a non- 
ideal experiment, it is important to remember that Equations (l34l l 
to (iJTl l come from the linear approximation, but Eq. ( l33l l does 
not. Hence, we expect our analytical calculations for the pres- 
sure to hold for much larger initial perturbations than the ones 
for the density. If the initial pressure disturbance is not small, 
but the linear expression for the plasma pressure is still valid, 
the adiabatic condition (i.e. conservation of entropy), which is 
more robust than the linear calculations, gives us a better ap- 
proximation for the final equilibrium plasma density, calculated 



as 



p(x,y, t — > oo) = 



p(x,y, f — » oo) p^(x,y, f = 0) 
p(x,y,t = 0) 



i/r 



(38) 



3. Numerical Experiments 

To investigate the validity of the analytical results, we have 
used Lare2D, a staggered Lagrangian-remap code with user con- 
trolled viscosity, to solve the full MHD equa tions, with the re - 
sistivity set to zero (for further details, see lArber et al.ll200ll) . 
The numerical domain is a square box with a uniform grid of 
256x256 points. The background magnetic field is pointing in 
the vertical y-direction and all the perturbations depend on both 
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X and y. The top and bottom boundaries of the box are periodic, 
so that the field lines are not line-tied. The boundaries on the left 
and right sides are closed. The choice of either closed or periodic 
side boundaries makes no difference for our experiments. There 
is neither mass nor energy flowing across the side boundaries. 

The numerical code uses the normalised MHD equations, 
where the normalised magnetic field, density and lengths, 

X — Lx , y — Ly , B = B„B , p — p„p , 

imply that the normalising constants for pressure, internal en- 
ergy and current density are 



y"0 



and ;„ = — - 



We have used here the subscript n for the normalising constants, 
instead of 0, to avoid confusion with the initial background 
values. The hat quantities are the dimensionless variables with 
which the code works. The expression for the plasma beta can 
be obtained from this normalisation as 



2p ^ 2pe(y - 1) 

£2 £2 



(39) 



1.5x10- 



1.0x10-" 



5.0x10- 




Total energy 
Internal energy 
Magnetic energy 
Kinetic energy 



1.5x10- 



1.0x10- 



5.0x10-^ 



hf^WvvvVWvvv^^ 




time 



Fig. 1. Time evolution of the energies for the case where 
max(/?i) = po = O.lB^/yLfQ andyS - 0.2, integrated over the whole 
2D box. The four evolutions have been shifted in the vertical axis 
by substracting the constant values 0, 0.5, 0.15 and 0.65, respec- 
tively, for kinetic, magnetic, internal and total energy, but their 
relative amplitudes are not scaled. The final losses of internal 
energy are entirely balanced with a net increase of magnetic en- 
ergy. The top figure is logarithmic in time and covers the whole 
relaxation. The bottom figure is linear in time and only covers 
the first part of the relaxation. From this graph, we can appre- 
ciate the complex oscillation periods, which are products of the 
sum of the diff'erent plane waves that drive the relaxation. 



The initial disturbance consists of an internal energy and 
pressure enhancement, leaving the initial plasma density and 
magnetic field undisturbed. The perturbation is taken as a single 
two-dimensional Gaussian centered in the middle of the box: 



ei(x, 0) = a exp 



(x - by 



exp 



(y-bf 

2c2 



(40) 



with b - 0.5L and c - 0.05L. As there is no initial perturbation 
of the magnetic field, the initial perturbed total pressure is just 
the perturbed plasma pressure, and the final analytical perturbed 
total pressure is the average value of that initial perturbed plasma 
pressure distribution, so that 

PT(oa) - pj(x) ^ J(^jpi(x,y,0)dy^ dx - J^pi(x,y,0)dy . 

For the experiment shown here, we have chosen a rather 
large perturbation, of the same order as the background value, 
i.e. a = eo, for which one may expect linear theory not to be 
applicable. The background plasma beta is 0.2. 

Figure [T] shows the time evolution of the various energies of 
the system, integrated over the whole box. Kinetic energy (in 
green) grows quickly from zero to its maximum value, and is 
subsequently damped to zero in the final equilibrium. A small 
fraction of the internal energy is converted into magnetic energy 
at the new equilibrium. For this particular set up, in which the 
perturbation has been introduced in the plasma pressure, it is the 
plasma that loses some of its initial internal energy, transfering it 
to the magnetic field. Note, that the amount of energy transfered 
is directly proportional to the magnitude of the perturbation. The 
total energy (i.e. the sum of the magnetic, kinetic and internal 
energies) is conserved to an accuracy of ~ 10"^. 

The 2D contour plots of the normalised plasma pressure, 
density and perpendicular current density for the initial and final 
equilibrium state are shown in Fig.|2]for this case. The initial in- 
crease of plasma pressure creates a localized decrease in plasma 
density. The displacement of the magnetic field lines is hard to 
appreciate from the field lines themselves, but is clear from the 
non-zero perpendicular electric current density, j^, that exists in 
the final equilibrium. 

Figures|3]and|4]show vertical cuts at x = L/2 and horizontal 
cuts aiy - Ljl, respectively, through the contour plots in Fig.|2] 
for the relevant quantities for the initial and final state. In these 
plots, a detailed comparison of the numerical and analytical so- 
lutions can be made. The analytical predictions for total pressure 
are very good, even though the amplitude of the perturbation is 
relatively large {a - eo) and so, strictly speaking, our analytical 
approximation should not hold. However, the linear approxima- 
tion for the plasma density (red crosses in Figures |3] and |4]i does 
not fit well. Instead, if we use Eq. ( l38T l, a much better fit for 
the plasma density is found (blue crosses in the plots), implying 
that the process is approximately adiabatic. Since the numerical 
experiments have been performed using a full MHD code that 
solves the non-linear equations, the process is not entirely adia- 
batic, but must have a finite amount of viscous heating that will 
become important as the initial perturbation is increased. 



4. Importance of non-linear effects 

To study how non-linearity affects the results as the magnitude 
of the initial perturbation increases, we focus again on the total 
pressure. The total pressure of the final numerical equilibrium 
must be constant, whether the relaxation remains in the linear 
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Fig. 2. Two-dimensional contour plots of plasma pressure (top), 
plasma density (center) and current density (bottom) in the initial 
state (left column) and final equilibrium (right column), for the 
same experiment as Fig. [1] 
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Fig. 3. Vertical cuts for the plasma pressure (left) and plasma 
density (right), for the same experiment as Fig. [1] and Fig. |2] 
Initial perturbed state (dashed) is compared with the final equi- 
librium, as found by the full MHD numerical simulations (solid) 
and predicted by the linear analysis (red crosses). For the den- 
sity predictions, the blue crosses represent predictions from the 
adiabatic condition given by Eq. ( |38] ). 



regime or not. On the other hand, the analytical definition of to- 
tal pressure given by ( fT3] l is an approximation from the linear 
analysis, and will become less valid as the non-linear terms be- 
come more important. We perform a series of experiments for 
various plasma beta values in which the relative amplitude of 
the initial perturbation is changed from a very small value, well 
within the linear regime, to a large value way outside it. Using 
these experiments, we investigate how the final total pressure de- 
parts from the linear predictions for different background plasma 
beta values. 

But first, we recall that the 2D relaxation may be separated 
into a vertical non-magnetic evolution (vertical redistribution of 
plasma pressure) and a horizontal evolution (horizontal redistri- 
bution of total pressure), in which the total pressure in the ver- 



0.10 
0.09 
= 0.08 
a 0.07 
0.06 
0.05 



0.0 0.2 0.4 0.6 

x/L 



1.0 



0.0 0.2 0.4 0.6 0.8 1.0 
x/L 



0.68 
CD^ 0.66 
S 0.64 
^ 0.62 

0.60 



0.0 0.2 0.4 0.6 O.E 
x/L 




0.0 0.2 0.4 0.6 0.8 1.0 
x/L 



Fig. 4. Horizontal cuts for the plasma pressure (top left), plasma 
density (top right), total pressure (bottom left) and magnetic field 
strength (bottom left), for the same experiment as Figures [T] |2] 
and[3] 

tical case is not determined by the linear analysis. This suggests 
that, effectively, in order to find a significant error in the final to- 
tal pressure for the 2D experiment, we will need very large val- 
ues of the initial two-dimensional perturbation. Hence, the fol- 
lowing experiments have been made for just a one-dimensional 
perturbation across the field lines. These results may be mapped 
onto those for our intital 2D perturbation, using the following: 



max(pi) 

PC) 



2D lexpi-iy-b)V2c^)dy 



max(pi)'^ 
Po I 



Figure |5] shows the relative error of the linear aproximation 
in both ID and 2D for the total pressure, as a function of the 
amplitude of the initial perturbation, for five different values of 
the plasma beta (J3 = 0.05, /? = 0.1, y6 = 0.2, /? = 1.3 andyS = 2). 
The bottom x-axis shows the magnitude of the one-dimensional 
perturbation, and the top .jc-axis shows the magnitude of the ini- 
tial two-dimensional perturbation before its vertical expansion. 
The error on the y-axis is calculated as the maximum difference 
between the linear prediction and the numerical results for the 
total pressure. 



max(|p'/" - p"/' 



1) 



where p™'" is the final constant total pressure obtained from the 
numerical simulations, and, for our non-magnetic initial pertur- 
bation. 



Po + Pl(x,oo) + - — + 

2/Uo juo 



with pi(x,oo) and By(x, oo) being the final perturbed plasma 
pressure and perturbed magnetic field from the numerical simu- 
lations. 

As j6 — » oo, we expect the relative error of the linear anal- 
ysis to tend to zero, independently of the perturbation, as in 
this case, the magnetic effects dissapear, and the initial pressure 
perturbation completely redistributes to a well defined constant 
value in the whole box. On the other hand, if /? <K 1, then the 
magnetic field will dominate over the plasma contributions, and 
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Fig. 5. Relative error in the linear prediction of the total pres- 
sure against the magnitude of the initial pressure perturbation, 
for five different values of the plasma beta. The slow growth rate 
of the error (non-linear eff'ects) indicates the validity of the linear 
analysis studied in this paper. 
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Fig. 6. Relative error of the density predicted assuming an adi- 
abatic evolution, with Eq. |38] against the magnitude of the 2D 
initial pressure perturbation, for three values of the plasma beta. 
Note, that the jc-axis in this plot is to be compared with the top 
.jc-axis in Fig.|5] 



much larger values for the relative initial perturbation will be 
needed for strictly leaving the linear regime. These two behav- 
iors can be seen in Fig.|5] where the plots for large plasma-betas 
tend to a smaller error, while the plots for small plasma-betas 
take longer to reach significant errors, i.e. to escape from the 
linear regime. Furthermore, we must not forget that we are here 
only talking about the initial background plasma beta, so a large 
background beta combined with a large initial perturbation will 
make the final plasma beta even higher. Thus, msx{p\)l po — > oo 
will imply — > oo for the final equilibrium, so we expect the 
curves of the relative error of the linear analysis to turn back to 
zero as the initial perturbation is greatly increased. In terms of 
energy conservation, as the velocity is zero at the initial and final 
states, the integral over the whole domain of internal energy plus 
magnetic energy must be conserved: If jS — > oo, then the internal 
energy is much larger than the magnetic energy, and will just re- 
distribute the plasma pressure, without transferring any energy 
into the magnetic field. 

On the contrary, the final plasma density is entirely deter- 
mined by the linear analysis, in both the vertical and the hori- 
zontal evolutions along and across the field lines, or in a better 
approximation, by the adiabatic condition. Hence, the non-linear 
effects for the plasma density will grow much quicker, as shown 
in Fig.|6l These last numerical experiments have been made for 
the original two-dimensional Gaussian perturbation. The error 
on the y-axis is given by 

maxOo'^'-p"""') 

where p"'' is the plasma density given by Eq. (l38T l, and p"'"" is 
the final density obtained with the numerical experiments. 

The relative error in the plasma density is considerably big- 
ger than the relative error in pressure, and so, for only a small 
change in pi / po in the 2D case, we find a large error in p. As this 
error quickly reaches significant values, the plasma beta plays 
much less of a role for the non-linear effects in the plasma den- 
sity than in the above total pressure. 



5. Summary and Conclusions 

We have presented analytical and numerical calculations for the 
2D magnetohydrodynamic relaxation of an untwisted perturbed 
magnetic system embedded in a plasma with beta of any size, 
and find a final equlibrium state which differs substantially from 
the initial background configuration. The equilibrium reached 
is a non force-free state in which the plasma pressure gradients 
are balanced by the magnetic Lorentz force. For a set of speci- 
fied boundaries, all the hydromagnetic quantities are fully deter- 
mined by the initial perturbed state. 

The initial disturbance evolves into the final relaxed state 
by different families of magnetoacoustic waves, dissipated via 
viscous damping. Fast magnetoacoustic waves propagate across 
the field lines in the horizontal direction, slow magnetoacoustic 
waves redistribute the thermodynamic quantities along the field 
in the vertical direction, and an extra contribution of slow mag- 
netoacoustic waves propagates along the magnetic field lines in- 
troducing a magnetic tension term. Nevertheless, these last slow 
magnetoacoustic waves dissipate the magnetic tension in such 
a way that it is totally unimportant when determining the fi- 
nal equilibrium distributions. The vertical redistribution of the 
plasma pressure to a homogeneous value demands the magnetic 
tension to dissapear completely, so both the plasma pressure and 
total pressure are one-dimensional at the end of the process. 

Within the linear regime, the final distributions are com- 
pletely independent of the viscosity, even though it is required 
to permit the relaxation to ocurr, as it is the only damping mech- 
anism of the waves. An increase in the viscosity enhances the 
diffusive term in the wave equation, and so, accelerates the pro- 
cess, but the final distribution is not modified. Instead, in the 
final equilibrium, all the quantities are simply determined by 
the behaviour of the final equilibrium total pressure, involving 
plasma and magnetic effects. Hence, the final equilibrium states 
for plasma pressure and magnetic field do not differ if the initial 
perturbation is of the density, temperature or internal energy. 

Finally, by investigating the linear regime, we have been able 
to make analytical predictions for the final MHS equilibrium, 
even when the regime is far from linear. The linear predictions 
remain remarkably valid even outside the linear regime, as the 
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growth rate of the non-ideal effects is very small, compared to 
the initial perturbations. 

In this paper, the introduction of plasma effects in the re- 
laxation of hydromagnetic systems have been studied in sim- 
ple schemes, producing a series of analytical predictions which 
have been confirmed by our numerical results. By starting with 
a uniform magnetic field, we have reached a state in which the 
magnetic field itself remains almost uniform, but where some 
current density has been built. Also, even in this simple con- 
figuration, a non-neglectible amount of energy has been trans- 
ferred from the plasma to the magnetic field. These implications 
of energy transfer during the relaxation process indicate that this 
process will be of importance in the solar corona, specially in 
the study of magnetic null points and their surroundings. Null 
points have been found to have a reas onable population density 
in the Solar Corona bv iLongcope & Parnell (2009). In further 
studies, we will consider more complex scenarios, introducing 
two-dimensional null points and starting to consider the impli- 
cations of non-zero plasma betas in three-dimensional magnetic 
environments. 
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